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Abstract 

We  compute  the  mutual  information  between  stored  and  retrieved  signals  in  a  spectral  hole-burning  optical  memory. 
Under  the  assumption  of  bitwise  data  detection  we  evaluate  the  maximum  achievable  storage  capacity  and  density  that  can 
be  obtained  for  both  time-domain  and  frequency-domain  architectures.  Noise  arising  from  crosstalk  in  both  the  time-  and 
frequency-domains  is  combined  with  the  material  shot  noise  to  arrive  at  optimal  system  designs  in  terms  of  material  volume 
and  number  of  bits  per  spatial  location.  Time-domain  results  demonstrate  the  existence  of  an  optimum  data  signal  strength 
and  the  analogous  frequency-domain  results  indicate  an  optimal  hole  depth.  We  find  that  realistic  information  storage 
densities  of  40  bits/fim3  are  possible  using  typical  material  parameters.  ©2000  Published  by  Elsevier  Science  B.V.  All 
rights  reserved. 

PACS:  42.79.Vb 

Keywords:  Optical  memory;  Optical  data  storage;  Spectral  hole  burning;  Information  theory;  Shot  noise;  Crosstalk  noise 


1.  Introduction 

Optical  storage  systems  play  an  important  role 
within  the  overall  memory  hierarchy.  Two-dimen¬ 
sional  (2D)  memories  based  on  optical  discs  (e.g., 
CDs  and  DVDs)  are  a  popular  medium  for  data 
distribution  and  archival  storage.  Three-dimensional 
(3D)  optical  storage  offers  larger  capacity  and  higher 
transfer  rates  as  compared  with  2D  approaches  and  is 
nearing  commercial  reality  with  many  recent  demon¬ 
strations  of  both  holographic  and  non-holographic 
systems  [1-4].  These  2D  and  3D  systems  realize 
their  advantages  by  efficient  use  of  spatial  degrees  of 
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freedom  exclusively.  Spectral-domain  memory  sys¬ 
tems,  in  addition,  use  wavelength  to  achieve  four-di¬ 
mensional  (4D)  data  storage.  These  4D  systems  may 
realize  storage  density  improvements  of  103  to  106 
times  that  of  3D  memories. 

Optical  storage  systems  based  on  spectral  hole- 
burning  (SHB)  materials  offer  the  high  data  rate  and 
low  latency  required  of  random  access  memory 
(RAM)  as  well  as  the  high  density  required  of  a 
large  mass  memory  (e.g.,  magnetic  disc).  Both  types 
of  system  have  been  studied  [5-7].  Both  pointwise 
(e.g.,  spatially  multiplexed)  and  volumetric  (e.g., 
holographic)  approaches  to  SHB  storage  have  been 
proposed  for  use  in  both  time-domain  and 
frequency-domain  architectures  [8,9].  Common  to  all 
of  these  SHB  systems  is  a  limit  on  the  smallest 
number  of  atoms  (or  molecules)  that  can  be  used  to 
store  a  single  bit  of  data.  As  storage  density  in- 
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creases  this  limit  is  approached,  and  material  shot 
noise  (i.e.,  uncertainty  concerning  the  actual  number 
of  atoms  participating  in  the  storage  process)  be¬ 
comes  important  [9,10].  Together  with  crosstalk 
among  stored  bits  (i.e.,  intersymbol  interference  in 
the  time-  and  frequency-domains),  this  noise  source 
sets  an  upper  limit  to  the  storage  density  of  SHB 
memory  systems. 

In  this  paper  we  treat  the  SHB  memory  system  as 
a  simple  communication  channel.  We  use  the  tools 
of  information  theory  to  estimate  the  storage  capac¬ 
ity  and  density  of  SHB  systems  in  the  presence  of 
material  shot  noise,  time/frequency-domain  crosstalk 
noise  and  envelope  effects.  The  mutual  information 
between  stored  and  retrieved  signals  is  measured 
under  the  assumption  of  bitwise  data  detection  in 
order  to  compute  the  information  storage  capacity  of 
the  SHB  memory.  We  study  this  capacity  as  a  func¬ 
tion  of  various  system  parameters  including  storage 
volume,  number  of  bits  per  location,  and  hole  depth. 
We  find  several  important  tradeoffs  underlying  the 
optimal  use  of  SHB  resources.  In  Section  2  we 
presen:  our  analysis  for  frequency-domain  storage 
and  in  Section  3  we  describe  the  time-domain  ap¬ 
proach. 

2.  Frequency-domain  storage 

The  frequency-domain  model  considered  here 
consists  of  a  SHB  material  that  exhibits  both  homo¬ 
geneous  and  inhomogeneous  broadening.  We  take 
the  homogeneous  and  inhomogeneous  linewidths  to 
be  2AyH  (Lorentzian)  and  2Avl  (Gaussian),  respec¬ 
tively.  A  frequency-agile  laser  is  used  to  saturate  the 
SHB  material  absorption  at  a  particular  spectral  loca¬ 
tion  vx,  thus  recording  a  single  bit.  Because  roughly 
M=  Av{/Aph  locations  are  spectrally  resolvable, 
we  approximate  the  storage  capacity  of  one  spatial 
location  as  C  =  M.  Fig.  1  depicts  the  absorption 
profile  of  one  such  spatial  location  within  a  SHB 
material  in  which  5  spectral  holes  have  been  burned 
to  record  the  9-bit  binary  word  101010101.  Retrieval 
of  this  spectrally  recorded  data  is  accomplished  by 
probing  the  SHB  material  at  each  of  the  9  data  bit 
frequencies  vf,  i  =  1, . . .  ,9.  Throughout  this  work  we 
will  limit  our  attention  to  the  storage  and  retrieval  of 
data  from  a  single  spatial  location. 


Fig.  1.  A  9-bit  binary-valued  word,  101010101,  is  stored  in  the 
Gaussian  inhomogeneously  broadened  profile.  The  T  bits  are 
represented  by  the  5  spectral  holes,  each  with  a  Lorenztian 
homogeneous  profile. 


We  first  describe  the  effect  of  material  shot  noise 
on  the  capacity  of  SHB  storage  [9].  Consider  a  single 
bit  to  be  stored  in  a  single  spectral  location.  As  the 
physical  volume  used  to  record  this  bit  decreases,  the 
number  of  atoms  that  can  participate  in  the  storage 
process  also  decreases.  A  Poisson  distribution  gov¬ 
erns  the  number  of  atoms  within  a  single  spectral 
channel.  In  the  absence  of  any  other  noise  sources, 
the  retrieved  ‘one’  and  ‘zero’  data  signals  are  there¬ 
fore  described  by  the  following  probability  density 


functions  (PDFs): 

j°o..(«)  s=e_Ao,~r 

(1) 

ajl  +  exp(2£aL) 

A°  2  Jd2Nx 

(2) 

A  01  • 

'  2'Jd2NA  exp(  £aL  -  1) 

(3) 

where  n  is  the  number  of  atoms  participating  in  the 
storage  in  a  single  channel,  <rs  is  the  absorption 
cross-section,  £  is  relative  hole  depth,  a  is  the 
absorption  efficiency,  Nk  is  the  density  of  the  atoms 
absorbing  wavelength  A  before  recording,  and  d  and 
L  are  the  spot  diameter  and  thickness  of  the  storage 
material,  respectively.  Fig.  2  depicts  example  PDFs 
for  two  different  values  of  the  storage  volume  with 
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Fig.  2.  Probability  density  functions  for  the  intensities  of  the 
retrieved  signals  0  and  1  for  two  different  values  of  the  storage 
volume:  (a)  V  =  2  \xm3\  and  (b)  V  =  20  |xm3.  (c)  Mutual  informa¬ 
tion  for  the  frequency-domain  storage  channel  characterized  by  a 
Poisson  intensity  probability  density  function.  The  mutual  infor¬ 
mation  is  a  function  of  spot  volume  for  single  bit  storage. 


a  =  10/cm,  £=  1,  Nx-  6.02  X  1013  mol/cm3  and 
A  =  600  nm.  Although  the  V  ==  20  |xm3  case  (Fig. 
2(b))  appears  well  separated,  we  observe  significant 
overlap  between  the  retrieved  1  and  0  signal  levels  in 
the  case  of  V  —  2  jxm3  (Fig.  2(a)).  (Note  that  re¬ 
trieved  data  bit  error  rate  is  proportional  to  the 
shaded  area  in  Fig.  2(a)  and  Fig.  2(b).)  This  result 
suggests  that  for  recording  volumes  approaching  the 
Rayleigh  volume  VR  =  A3,  the  retrieved  signal  fi¬ 
delity  can  be  quite  poor  owing  to  material  shot  noise. 

It  is  possible  to  quantify  the  capacity  impact  that 
arises  from  the  loss  in  retrieved  signal  fidelity  de¬ 
scribed  above.  When  data  is  retrieved  from  a  noisy 
memory  system,  each  stored  bit  gives  rise  to  less 
than  one  retrieved  bit.  The  uncertainty  that  arises 
from  noise  manifests  itself  as  a  loss  of  information. 
This  loss  of  information  is  characterized  by  the 


mutual  information  R  between  the  stored  and  re¬ 
trieved  signals: 


R  =  max 

(TT0,7r,) 


x  log  2 


PoiO 

V oPo( 7)  +  P\(J) 


d  / 


+  IT 


./  />■(/) 

J  —  3C 


Xlog2 


PjJ) 

KoPoi1)  +  ^iPli1) 


(4) 


where  p0(I)  and  p  ,( /)  are  PDFs  for  the  retrieved 
signals  0  and  1,  and  7r0  and  vl  are  the  prior 
probabilities  for  0  and  1  data  bits,  respectively  [11]. 
The  mutual  information  R  therefore  represents  the 
actual  information  capacity  of  a  single  spectral  chan¬ 
nel  of  the  memory  system.  Fig.  2(c)  shows  a  plot  of 
mutual  information  versus  volume  for  the  single  bit 
SHB  storage  system  described  above.  As  material 
volume  decreases,  the  information  per  bit  decreases 
reflecting  the  signal  degradation  arising  from  shot 
noise;  while,  for  large  volumes  the  mutual  informa¬ 
tion  approaches  1  indicating  that  noise  has  become 
negligible. 

Now  we  would  like  to  extend  this  analysis  to  the 
case  of  multiple  spectral  locations.  We  consider 
recording  M  bits  within  the  inhomogeneously 
broadened  SHB  material  profile.  In  this  case  a  frac¬ 
tion  A  vx/M  of  the  spectrum  can  be  allocated  to 
each  bit.  Given  a  fixed  material  volume  it  is  once 
again  possible  to  compute  the  number  of  active 
atoms  participating  in  the  storage  process  for  each 
bit  and  thus  the  fidelity  of  memory  recall.  We  expect 
a  tradeoff  governed  by  the  finite  SHB  material  re¬ 
sources.  As  M  increases,  the  number  of  active  atoms 
participating  in  the  storage  process  for  each  bit  must 
decrease.  This  results  in  a  storage  capacity  tradeoff 
related  to  balancing  an  increase  in  the  number  of  bits 
against  the  concomitant  loss  of  information  per  bit. 
We  will  once  again  use  the  mutual  information  to 
quantify  this  tradeoff  as  C  =  MR(M\  where  C  is  the 
storage  capacity,  and  the  mutual  information  R(M\ 
has  been  made  an  explicit  function  of  the  number  of 
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stored  bits.  Shot  noise  however,  is  not  the  only  error 
source  that  affects  R(M).  Three  additional  effects 
must  also  be  discussed. 

The  first  additional  noise  source  to  be  included  is 
spectral  crosstalk  among  nearby  bit  locations.  Be¬ 
cause  the  spectral  resolution  is  governed  by  the  SHB 
material  homogeneous  broadening,  the  tails  of  the 
corresponding  Lorentzian  distributions  will  result  in 
crosstalk  among  the  stored  data.  As  the  number  of 
spectral  location  is  increased,  this  crosstalk  noise 
becomes  more  severe  representing  another  contribu¬ 
tion  to  the  tradeoff  between  M  and  R(M).  In  our 
model  the  areas  under  the  tails  of  each  4  nearest 
neighbors  Lorentzian  distribution  function  are  used 
to  estimate  the  crosstalk  noise.  The  crosstalk  at  the 
ith  channel  with  central  frequency  vt  is  determined 
by 


i  +  2 

x(i)=  £  xu 

j-i- 2 
j*< 


(5) 


v  fvi  +  *VJ  I  I  \A  /••',  +  A*' ,/M 

x>r)  .  dMv)dv=J  .  /udj 

’'pj  —  Av  J  Avx/M 


*(")- 


(a^h/2  y 


( v-Vj )  +(Avh/2) 


g(v) 


du 

(6) 


where  Xu  is  the  crosstalk  at  the  spectral  channel  vt 
due  to  the  channel  at  vj3  g{v)  is  the  Gaussian 
inhomogeneous  profile,  Lj(v)  is  the  Lorentzian  ho¬ 
mogeneous  profile  at  and  d-  is  the  data  0  or  1 
recorded  at  channel  j.  The  second  additional  noise 
effect  is  also  related  to  crosstalk  noise  and  concerns 
the  impact  that  hole  depth  will  have  on  retrieved 
signal  fidelity.  All  of  the  results  discussed  to  this 
point  have  utilized  hole  depth  £=  1,  representing  the 
largest  possible  retrieved  signal  level  and  therefore 
the  smallest  possible  shot  noise.  Large  £  however, 
can  exacerbate  crosstalk.  Because  small  hole  depth 
corresponds  to  large  shot  noise  while  large  hole 
depth  produces  excess  crosstalk,  we  expect  to  see  a 
tradeoff  between  £  and  C  resulting  in  an  optimum 
value  of  hole  depth.  The  relation  between  retrieved 


signal  strength  and  a  variety  of  material  parameters 
including  hole  depth  is  given  by: 

/(£)  =  -^-[exp(-aL(l  -  f  ))  -exp(-aL)] 

(7) 

where  q  is  the  elementary  charge,  rjQ  is  the  detector 
quantum  efficiency,  P  is  the  readout  laser  power, 
and  hv  is  the  photon  energy  [9].  The  third  effect  to 
be  included  in  our  model  is  nonuniform  fidelity  of 
holes  as  a  function  of  spectral  location.  From  Fig.  1 
it  is  clear  that  holes  recorded  near  the  SHB  absorp¬ 
tion  peak  at  v0  utilize  a  larger  number  of  active 
atoms  than  those  recorded  near  the  edge  of  the 
inhomogeneous  spectrum.  This  nonuniformity  will 
result  in  additional  spreading  of  the  retrieved  signal 
PDFs  and  a  further  decrease  in  the  mutual  informa¬ 
tion.  We  have  used  a  Gaussian  model  to  represent 
the  inhomogeneous  profile  and  to  compute  the  result¬ 
ing  retrieved  signal  PDFs. 

Combining  the  three  noise  sources  described 
above  with  the  effect  of  material  shot  noise  de¬ 
scribed  earlier,  we  compute  the  SHB  storage  capac¬ 
ity  as  a  function  of  M  for  a  fixed  storage  volume 
and  several  values  of  £.  The  results  of  these  calcula¬ 
tions  are  shown  in  Fig.  3.  From  this  figure  we  see 
that  for  a  storage  volume  of  V  -  8  |xm3  the  storage 


Fig.  3.  Capacity  versus  the  number  of  spectral  channels  for 
various  values  of  hole  depth  f  at  a  given  spot  volume  V  —  8  p,m3. 
Both  large  (£  —  0.8)  and  small  (f  =  0.2)  hole  depths  result  in  a 
lower  capacity  than  the  moderate  hole  depth  (  £  =  0.6). 
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capacity  initially  increases  rapidly  as  the  number  of 
spectral  locations  is  increased.  Eventually  however, 
shot  noise  and  crosstalk  begin  to  dominate  so  that 
further  increases  in  M  are  offset  by  losses  in  R(M). 
This  effect  gives  rise  to  slowly  decreasing  capacity 
behavior  observed  beyond  some  optimum  number  of 
spectral  channels.  We  also  see  that  both  small  (  £  — 
0.2)  and  large  (£  =  0.8)  values  of  hole  depth  result 
in  lower  capacity  than  the  moderate  value  of  £  =  0.6. 
This  is  consistent  with  the  expected  tradeoff  involv¬ 
ing  hole  depth.  Extracting  the  capacity  peaks  from 
Fig.  3  and  plotting  these  optimal  capacity  values  as  a 
function  of  hole  depth,  we  obtain  the  curve  labeled 
V  =  8  pirn3  shown  in  Fig.  4.  The  existence  of  an 
optimum  hole  depth  is  obvious  from  this  figure.  Fig. 
4  also  presents  the  optimum  capacities  for  two  other 
values  of  storage  volume  and  we  see  that  capacity 
increases  with  storage  volume.  This  is  to  be  expected 
owing  to  the  additional  storage  resources  that  come 
from  increasing  the  storage  volume.  An  important 
question  however  concerns  how  best  to  use  this 
volume.  Figs.  3  and  4  facilitate  optimal  system  de¬ 
sign  with  respect  to  both  the  number  of  spectral 
locations  and  the  hole  depth.  Extracting  these  opti¬ 
mum  designs  from  Fig.  4  and  dividing  by  the  associ¬ 
ated  storage  volumes  results  in  the  storage  density 
curve  shown  in  Fig.  5.  Although  capacity  increases 
with  volume,  high  density  requires  both  high  capac¬ 
ity  and  small  volume,  resulting  in  the  tradeoff  de¬ 
picted  in  this  figure.  An  optimum  storage  density  of 
D0 pt  =  40  bits/jxm3  is  obtained  at  an  optimum  vol- 


Fig.  4,  Optimum  capacity  versus  hole  depth  f  for  three  different 
spot  volumes:  V  =  8  p.m3,  V  =  24  |xm\  and  V  =  32  p,m3. 


Fig.  5.  Information  density  is  a  function  of  the  spot  volume.  We 
obtain  the  optimum  density  of  40  bits/fxm3  at  a  volume  VopI  =  9 
fim3,  a  hole  depth  of  fopt  —  0.6,  and  the  number  of  spectral 
channels  Mopl  =  650. 

ume  of  Vopl  =  9  p.m3,  an  optimum  hole  depth  of 
f  =  0.6,  and  the  optimum  number  of  spectral  chan¬ 
nels  Mopt  =  650. 

These  results  suggest  that  previous  estimates  of 
SHB  storage  densities  were  overly  optimistic.  They 
were  optimistic  with  respect  to  both  the  number  of 
spectral  channels  and  the  storage  volume.  In  particu¬ 
lar,  we  have  assumed  a  ratio  Av{/AvH  =  1000  yet 
the  optimum  number  of  spectral  channels  in  the 
presence  of  shot  and  crosstalk  noise  is  found  to  be 
Afopt  =  650.  Similarly,  the  Rayleigh  volume  for  our 
model  is  VR  =  0.216  |xm3;  however,  the  volume  that 
optimizes  storage  density  was  found  to  be  Vopt  =  9 
fim3,  nearly  42  times  the  Rayleigh  volume. 

Application  of  the  frequency-domain  approach  to 
SHB  storage  has  been  limited  due  to  the  potentially 
severe  latency  costs  associated  with  high  spectral 
resolution.  For  example,  the  spectral  features  corre¬ 
sponding  to  the  optimal  system  described  above 
result  in  single  bit  latencies  in  the  range  of  r  = 
1  /A  vH  =  1  (jls.  Readout  of  a  single  650-bit  location 
using  a  single  frequency  tunable  laser  requires  a  time 
of  650  pus  and  is  therefore  unacceptable  for  most 
applications.  An  alternative  to  the  frequency-domain 
approach  is  the  use  of  pulsed  lasers  in  a  time-domain 
storage  architecture  based  on  SHB.  Such  a  system 
solves  the  latency  problem  by  reading  all  spectral 
channels  simultaneously;  however,  similar  noise  is¬ 
sues  arise  in  the  time-domain  approach.  We  investi¬ 
gate  these  issues  in  the  next  section. 
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3.  Time-domain  storage 

We  can  consider  time-domain  storage  to  be  an 
application  of  spectral  holography,  storing  each  bit 
not  in  any  one  spectral  location  but  instead  on  all 
spectral  locations.  Reference  and  data  laser  pulses 
are  mixed  to  produce  a  temporal  interference  pattern, 
which  is  recorded  on  the  SHB  material  population. 
To  illustrate  this  process,  we  refer  to  Fig.  6  to 
describe  the  storage  and  retrieval  of  a  data  train.  Fig. 
6(a)  and  Fig.  6(b)  show  the  storage  and  retrieval 
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Data  pulse  train 
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Fig.  6.  A  typical  SHB  storage  system  using  the  time-domain 
architecture:  (a)  the  interference  pattern  between  reference  and 
data  pulses  is  stored  in  the  SHB  material  during  the  recording 
process;  (b)  a  data  echo  modulated  by  the  exponential  dephasing 
is  retrieved  after  a  read  pulse  passes  through  the  material;  and  (c) 
the  interference  pattern  is  recorded  by  modifying  the  ground-state 
population  of  the  material. 


processes,  and  Fig.  6(c)  depicts  the  ground  state 
population  of  the  SHB  material  upon  reacting  to  the 
interference  of  two  temporal  sequences.  The  refer¬ 
ence  pulse  should  be  sufficiently  brief  to  uniformly 
excite  the  storage  material  throughout  its  entire  inho¬ 
mogeneous  bandwidth.  As  the  reference  pulse  r(t) 
passes  through  the  material,  a  coherence  in  the  stor¬ 
age  material  represented  by  R(v)  is  generated,  where 
R(v)  is  the  temporal  Fourier  transform  of  r(r).  If  the 
data  pulse  train  d(t)  interacts  with  the  material 
before  this  coherence  has  relaxed,  the  material  reacts 
to  the  superposition  of  these  two  spectra  by  modify¬ 
ing  the  ground-state  atomic  population  of  the  mate¬ 
rial.  For  recording  in  the  linear  regime,  the  induced 
ground  state  population  is  approximately  given  by: 

Pg(v)  a\R(v)  +  D( j')exp(i2iri'0T)|2,  (8) 

where  Pg(v)  is  the  population  of  the  SHB  atoms  in 
their  ground  state,  D(v)  is  the  Fourier  transform  of 
the  data  signal  d(t\  and  r  is  the  delay  between  the 
reference  and  data  signals.  This  interference  between 
reference  and  data  signals  is  a  spectral  hologram  that 
can  be  retrieved  by  using  a  read  pulse  having  the 
same  temporal  form  as  the  reference  pulse.  Suppose 
the  material  is  illuminated  by  the  read  pulse  at  some 
later  time.  This  readout  process  will  induce  coherent 
radiation,  known  as  a  stimulated  photon  echo  and  a 
train  of  signal  pulses,  es(t\  is  obtained  which  has 
roughly  the  same  temporal  profile  as  the  original 
data  pulse  train, 

es(t)  a  f  fLexp(i2Trvt)R*  (v)D(v) 

1  —cc  •'O 

XR(  v)g{  v)exp[  —  a(  v)l]dldv.  (9) 

It  should  be  pointed  out  that  Eq.  (9)  is  only  valid  for 
aL  «:  1  and  indeed,  this  condition  is  satisfied  for  all 
of  the  results  presented  in  this  paper. 

An  estimate  of  the  storage  capacity  of  the  time- 
domain  memory  can  be  obtained  by  considering  the 
constraints  placed  on  the  data  signal  by  the  proper¬ 
ties  of  the  SHB  material.  The  first  constraint  is  that 
the  data  pulse  bandwidth  should  not  exceed  the  SHB 
material  bandwidth,  which  suggests  the  requirement 
8f>  1/Avj.  The  second  restriction  concerns  the 
duration  of  the  data  train,  which  is  determined  by  the 
dephasing  effect  according  to  the  Lorentzian  spectral 
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profile  in  the  frequency-domain.  This  constraint  re¬ 
sults  in  a  maximum  time  between  the  first  and  last 
data  pulses  given  by  T <  l/kvH.  Thus,  approxi¬ 
mately  T/ht  —  Ai/j/A  vH  bits  of  information  can  be 
stored  in  a  single  spatial  location,  which  is  as  the 
same  as  the  potential  capacity  of  the  frequency-do¬ 
main  storage. 

Once  again  we  consider  the  effect  of  material  shot 
noise  and  begin  with  the  storage  and  retrieval  of  a 
single  data  bit.  When  the  material  is  exposed  to  brief 
reference  and  data  pulses,  their  relative  amplitudes 
determine  the  modulation  depth  of  the  spectral  inter¬ 
ference  pattern,  defined  as 


m(v) 


D(v)R*(v)+D’(v)R(v) 

|DO)l2  +  l*0)l2 


(10) 


where  R*(v)  and  D*(v)  are  the  conjugate  Fourier 
transforms  of  the  reference  pulse  and  data  pulse 
train,  respectively.  Frequencies  with  small  m(v )  suf¬ 
fer  more  severe  shot  noise  because  fewer  active 
atoms  are  involved  in  the  storage  process.  On  the 
other  hand,  some  frequencies  will  saturate  their  mod¬ 
ulation  depth  and  will  produce  distortion  in  the 
retrieved  data  signal  train.  This  implies  that  there 
exists  a  tradeoff  between  saturation  distortion  and 
shot  noise,  and  this  tradeoff  can  be  accessed  through 
the  modulation  depth  m(  v).  To  simplify  our  simula¬ 
tions,  the  reference  amplitude  is  fixed  to  satisfy  the 
requirement  of  optimized  tt/2  pulse  area  [6,12]. 
Then  m(v)  is  totally  determined  by  the  data  pulse 


Fig.  8.  Capacity  versus  the  number  of  data  bits  for  various  values 
of  data  pulse  amplitude  at  a  given  spot  volume  V  =  32  jim3.  Both 
large  (140  V/m)  and  small  (60  V/m)  £d  result  in  a  lower 
capacity  than  the  moderate  £d  (80  V/m). 

amplitude  £d.  Similar  to  the  analysis  for  frequency- 
domain  storage,  we  use  a  Poisson  distribution  to 
describe  the  number  of  excited-state  atoms  in  each 
spectral  channel  of  the  recorded  modulation  depth 
function.  We  compute  the  mutual  information  be¬ 
tween  the  stored  and  retrieved  signals,  and  plot  the 
dependence  of  mutual  information  on  data  pulse 
amplitude  in  Fig.  7.  As  the  data  amplitude  decreases, 
the  information  per  bit  decreases  reflecting  the  loss 
of  signal  fidelity  arising  from  increased  shot  noise. 
While  for  large  m(v)  the  mutual  information  also 
decreases,  owing  to  the  effect  of  saturation  distor- 
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Fig.  7.  Mutual  information  as  a  function  of  data  pulse  amplitude. 
There  exists  an  optimum  £d  balancing  the  material  shot  noise 
against  the  saturation  distortion. 


Fig.  9.  Optimum  capacity  versus  data  pulse  amplitude  for  three 
different  spot  volumes:  V  =  8  fxm3,  V  =  24  fim3,  and  V  —  32 
p,m3. 
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Fig.  10.  Information  density  versus  spot  volume.  We  obtain  the 
optimum  density  of  42  bits/|im3  at  a  volume  Vopt  =  8  p,m\  a 
data  pulse  amplitude  of  £d  opl  =  80  V /m,  and  the  number  of 
spectral  channels  Mop{  =  750. 

tion.  There  exists  an  optimum  data  amplitude  that 
balances  these  two  error  sources.  This  effect  is  anal¬ 
ogous  to  the  existence  of  an  optimum  hole  depth  for 
the  frequency  domain  case. 

This  analysis  can  be  extended  to  the  case  of 
storing  multiple  bits.  Dephasing  results  in  an  expo¬ 
nents 1  temporal  envelope  modulating  the  retrieved 
data  signals  [8].  In  addition,  temporal  crosstalk  among 
neighboring  bits  is  produced  by  the  Gaussian  inho¬ 
mogeneous  profile.  We  treat  the  storage  system  as  a 
simple  communication  channel  with  a  Gaussian  im¬ 
pulse  response.  Thus  the  output  signal  sequence  is  an 
exponential  modulation  of  the  convolution  of  the 
input  data  train  with  a  Gaussian  impulse  response. 
This  implies  that 

e,(t)  =  [d(t)®  G(t)]exp(-t/T2)  (11) 

where  G(t )  is  the  inverse  Fourier  transform  of  the 
inhomogeneous  profile  g(v).  In  the  presence  of  the 
three  effects  described  above:  material  shot  noise, 
crosstalk  noise,  and  the  dephasing  envelope,  we 
simulate  the  storage  and  retrieval  process  again  and 


compute  the  mutual  information  of  each  bit  R(M). 
The  storage  capacity  is  determined  by  C  —  MR(M). 
Once  again  as  M  increases,  crosstalk  and  shot  noise 
increase,  yielding  a  tradeoff  between  the  increasing 
number  of  bits  and  the  decreasing  information  per 
bit.  Fig.  8  shows  capacity  C  versus  the  number  of 
data  bits  M  for  various  values  of  the  data  sequence 
amplitude  Ed  subject  to  a  fixed  storage  volume. 
From  this  figure  we  can  see  that  for  the  storage 
volume  of  V  =  32  fim3,  the  storage  capacity  initially 
increases  as  the  number  of  data  bits  increases.  But 
for  large  M,  the  crosstalk  and  shot  noise  eventually 
dominate  and  losses  in  R(M)  offset  the  increase  in 
M.  It  is  also  observed  that  both  small  ( Ed  —  60 
V/m)  and  large  (Ed  —  140  V/m)  data  amplitude 
result  in  lower  capacity  than  the  moderate  amplitude 
(£d  =  80  V/m).  This  result  is  consistent  with  the 
anticipated  tradeoff  involving  data  amplitude  £d,  and 
also  agrees  with  the  hole  depth  effect  seen  in  the 
frequency-domain  case.  Extracting  the  capacity  peaks 
from  Fig.  8,  we  plot  optimum  capacity  versus  data 
pulse  amplitude  in  Fig.  9  for  various  active  volumes. 
It  is  clear  to  see  an  optimal  data  amplitude  for  each 
storage  volume.  Extracting  the  optimal  capacity  val¬ 
ues  from  Fig.  9  and  dividing  by  the  corresponding 
storage  volumes  results  in  the  density  curve  shown 
in  Fig.  10.  An  optimal  storage  density  of  D  —  42 
bits/jjim3  is  obtained  at  an  optimal  data  amplitude 
of  80  V /m  in  an  optimum  volume  of  Vopt  =  8  jxm3, 
allowing  Mopt  —  750  bits.  Obviously,  this  density 
tradeoff  is  analogous  to  the  frequency  domain  case 
and  the  results  are  similar. 

Comparing  the  simulation  results  for  the  time  and 
frequency  domain  memories,  there  are  no  significant 
differences  in  storage  capacity  between  these  two 
architectures.  To  compare  the  frequency-  and  time- 
domain  storage  systems,  we  list  the  physical  parame¬ 
ters  in  Table  1.  But  the  intrinsic  advantage  of  the 
time  domain  memory  over  frequency  domain  mem- 


Table  1 

Parameters  used  for  simulations 


Frequency-domain 

Time-domain 

Physical  material 

A  vH  homogeneous  linewidth,  1  MHz 

Avl  inhomogeneous  linewidth,  1  GHz 
f  hole  depth 

M  number  of  spectral  channels 

8 1  duration  time  of  a  single  bit,  50  ns 

T1  duration  time  of  data  train,  50  |xs 
m(v)  modulation  depth 

M  number  of  data  bits 

C  concentration,  10 ~4  mol/1 

apeak  absorption  effeciency,  lOmol/cm 

Nx  No.  atoms  absorbing,  A/V 
d  spot  diameter,  2  \xm 
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ory  is  its  ability  to  store  and  retrieve  spectral  infor¬ 
mation  in  parallel.  Time  domain  storage  offers  a 
means  of  achieving  high  storage  density  and  high 
data  rate  simultaneously.  In  addition,  it  eliminates 
the  necessary  requirement  of  a  very  precise  fre¬ 
quency-agile  laser  (or  an  array  of  such  lasers). 

4.  Conclusions 

We  have  described  a  formalism  within  which 
information  theory  can  be  used  to  estimate  the  stor¬ 
age  capacity  and  volumetric  storage  density  of  spec¬ 
tral  hole-burning  memory.  We  have  used  this  formal¬ 
ism  to  analyze  both  time-  and  frequency-domain 
architectures.  Our  simplified  analyses  have  included 
the  deleterious  effects  of  material  shot  noise,  crosstalk 
noise,  and  both  time-domain  (exponential)  and  fre¬ 
quency-domain  (Gaussian)  envelope  functions.  We 
characterize  the  fidelity  of  data  retrieved  from  SHB 
memory  in  terms  of  the  information  per  bit.  Our 
results  reflect  an  important  tradeoff  between  the 
number  of  bits  stored  in  each  spatial  location  and  the 
retrieved  information  per  bit.  We  analyze  this  trade¬ 
off  in  both  the  time-  and  frequency-domains  to  find 
optimal  values  of  hole  depth,  data  pulse  amplitude, 
storage  volume,  and  the  number  of  spectral  channels. 
An  optimum  storage  density  of  40  bits/|xm3  is 
consistent  for  both  time-  and  frequency-domain  ar¬ 
chitectures. 

There  are  three  important  assumptions  that  under¬ 
lie  this  work.  The  first  concerns  our  focus  on  bitwise 
data  retrieval.  This  assumption  precludes  the  detec¬ 
tion  of  SHB  output  signals  using  advanced  signal 
processing,  equalization,  or  other  soft-decision  meth¬ 
ods.  Ongoing  work  is  focused  on  developing  new 
|  interface  processing  algorithms  for  use  in  SHB  mem¬ 
ories  and  the  inclusion  of  such  techniques  will  im- 
|  prove  the  capacity  and  density  estimates  reported 
herein.  Our  second  assumption  concerns  the  rela¬ 
tively  low  doping  density  (Nx  =  10,3/cm3)  used 
;  throughout  this  work.  This  assumption  is  relevant  to 
organic  systems  and  results  in  material  shot  noise 
playing  an  important  role  in  determining  the  re¬ 
trieved  signal  fidelity.  In  order  to  estimate  the  range 
of  concentration  for  which  our  assumption  is  valid, 
we  can  compare  the  material  shot  noise  with  thermal 
detector  noise.  Consider  a  candidate  high  speed  de¬ 


tector  and  amplifier  electronics  with  noise  equivalent 
power  equal  to  1  X  10" 10  W.  If  such  a  detector 
system  operates  at  a  data  rate  of  1  X  108  bits/s,  then 
the  noise  associated  with  the  detector  will  become 
larger  than  material  shot  noise  for  concentrations 
greater  than  1.2  X  10“3  mol/1  ( Nx  >  5  X  1014/ 
cm3).  For  such  heavily  doped  systems,  the  detector 
noise  must  be  included  to  obtain  accurate  informa¬ 
tion  theoretic  density  estimation  owing  to  the  rela¬ 
tively  low  oscillator  strength  associated  with  these 
inorganic  systems.  The  third  assumption  concerns 
our  focus  on  a  single  spatial  location  and  can  be 
relaxed  by  including  the  effects  of  beam  nonuniform¬ 
ities  and  spatial  crosstalk  within  the  computation  of 
the  mutual  information  per  bit. 
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